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Abstract 

We consider nonlinear parabolic SPDEs of the form dtu — — (— A) a//2 tt+ 
b(u) + a(u)w, where w denotes space-time white noise. The functions b 
and a are both locally Lipschitz continuous. Under some suitable condi- 
tions on the parameters of these SPDEs, we show that the second moment 
of their solutions blow up in finite time. This complements recent works 
of Khoshnevisan and his coauthors; see for instance [7], [8] and [5] as well 
as those of Chow [3] and [4]. Furthermore, upon comparing our stochastic 
equations with their deterministic counterparts, we find that our results 
indicates that the presence of noise might affect the occurrence of blow-up. 
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1 Introduction 

In [7], [8] and [5], Khoshnevisan and his coauthors initiated a research program, 
where they systematicaly studied the large time behavior of stochastic equations 
of the following type 

dtu(t, x) — £u(t, x) + a(u(t, x))F(t, x) ^ 
u(0 ,x) = u (x), 

where £ is the generator of a Levy process, a is a globally Lipschitz function, 
F is a noisy perturbation, and uq(x) is specified initial data satisfying certain 
suitable conditions. Among other things, the authors studied the upper p-th 
moment Liapounov exponent 7(p) defined as 

j(p) := limsup - lnE(|«(t, x)\ p ) for all p G (0, oo). 

t— >oo t 
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A motivation for this study comes from the Parabolic Anderson model and 
the phenomenon of intermittency which this model exhibits. Instead of going 
into details about this phenomenon, we refer the reader to [2] for details and 
references. One of the main results in the above works, is the existence and 
strict positivity of those upper Liapounov exponents. A common assumption in 
all the works mentioned above, is that the operator C models regular diffusion 
so one considers say the Laplace operator. Another important assumption is 
that the multiplicative non-linearity a is assumed globally Lipschitz. Indeed, 
this is required to ensure global existence of solutions. In this short note, we 
consider a subclass of (1.1), where we consider a fractional Laplacian operator 
instead of the diffusion process, and we drop the global Lipschitz property of a. 
More precisely we look at the following equation: 

d t u(t, x) = -(-A) a / 2 u(t, x) + b(u(t, x)) + a(u(t, x))w{t, x) 
u(0 ,x) = u (x), 

where {w(t, i)}(>o )ie n denotes space-time white noise. -(-A)"/ 2 is the frac- 
tional Laplacian; the L 2 -generator of a symmetric stable process of order a, so 
that Eexp(i£ • X t ) — cxp(— t|£| Q ). We will follow Walsh [18] and interpret the 
above as an Ito-type stochastic PDE. The functions b and a are assumed to 
be positive and locally Lipschitz. A standard localisation argument show that 
the above equation has a unique random field solution with finite moments, 
provided that 



confer with [7] and [6]. This implies that a e (1,2], a condition which will be 
enforced throughout this paper. We will also assume the initial data uq is a 
positive bounded measurable function satisfying the following lower bound: 



where k is a positive constant. 

There are many physical motivations to consider a fractional Laplacian op- 
erator. These operators appear in many models in non-Newtonian fluids, in 
models of viscoelasticity such as Kelvin- Voigt models, various heat transfer pro- 
cesses in fractal and disordered media and models of fluid flow and acoustic 
propogation in porous media. Interestingly, they have also recently been ap- 
plied to pricing derivative securities in financial markets. Sec [12, 11, 1] for 
details on a number of the afforcmentioned applications. The essential feature 
behind many of these models, is that they aim to capture the nonlinear relation 
between stress and strain, in non-Newtonian or viscoelastic models perse. Or 
perhaps, aim to model the disparity between scaling laws associated with dif- 
fusion, in cuclidian geometries, and fractal geometries, such as certain porous 




for all p > 0; 



inf u (x) > k, 



(1.3) 
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media. Thus there are a plethora of reasons why one might consider investigat- 
ing the asymptotic behavior of equations of the type (1.2). In particular, such 
investigations might help us glean information in many real physical scenarios, 
involving sub or super diffusions. 

The primary contribution of the current manuscript is to show that under 
certain conditions on the initial data and parameters, the second moment of the 
solution to (1.2), blows up in finite time, and hence we have non-existence of the 
upper Liapounov exponents. When we say that a function / blows up in finite 
time, we mean that there exists a time T such that for all t >T, f(t) = oo. And 
when we say that a statement holds whenever a quantity Q is large enough, we 
mean that the statement holds whenever Q > Qq, where Qq is some positive 
constant. 

Our result will hold under at least one of the following two conditions which 
are non-linear growth conditions on a and b respectively. 

Condition 1.1. There exist constants j5 and K\ > 0, both strictly positive such 
that 

inf ^ > K X . 

xeR \x\ 1+ p 

Condition 1.2. There exist constants 7 and K2 both strictly positive such that 

xeR \x\ 1+ i ~ 

Denote {u(t, i)}t>o,ieR as the solution to (1.2) so that {E|u(i, x)\ 2 } t >o, x eR 
denotes the corresponding second moment. 

Theorem 1.3. Suppose {u(t, x)} t >o. x eR denote the unique solution to (1.2). 
We then have the following 

(a) Under Condition 1.1, {E\u(t, 2;)| 2 }t>o, xeR the second moment of the so- 
lution, blows up in finite time whenever k is large enough. 

(b) Under Condition 1.2, {E\u(t, x)\ 2 } t >o, xeR blows up infinite time for any 
initial data satisfying (1.3), whenever 7 < a. But if 7 > a, finite time 
blow up occurs only whenever k is large enough. 

We remark that blow-up occurs pointwise, that is for each x £ R. This is in 
contrast with the results in [3] and [4] . 

Blow-up problems in deterministic cases, is an area of intense study, starting 
with the classical work of Fujita [9]. Blow up can mean the solutions them- 
selves tends to infinity in finite time or sometimes their derivatives blow up 
in finite time; for instance solutions to the Burgers' equation, starting with 
smooth initial data, can experience loss of smoothness in finite time [19], that is 
lim t ^T*<oo ||Vit||2 = +00. The review articles of [13] and [14] contain a detailed 
exposition of the subject. A common approach in tackling these various blow 
up problems is to consider a "right" quantity, which blows up in finite time, 
usually by exercising a differential inequality. The quantity at hand is usually 
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linked to a certain "norm" of the solution. We will use a similar approach to 
prove our main result. One of the main difficulties here is to find this "right" 
quantity; see Proposition 3.1. 

We now briefly give an outline of the paper. In Section 2, we recall some 
facts about the heat kernel of the fractional Laplacian and prove some important 
estimates. In Section 3, we prove the main result. Finally in section 4, we 
provide some concluding remarks. Throughout the paper, we use the letter c 
with or without subscripts to denote a constant whose value is not important 
and may vary from places to places. And by f(x) x g{x), we mean that there 
exists a constant K such that -^g(x) < f(x) < Kg(x). 



2 Some auxiliary results. 

We begin this section by developing some required background. Let p(t, •) 
denote the transition density function of a symmetric stable process of order a. 
The existence of this density is a well known fact; see for instance [15]. Set 



V t f(x):= / p(t,x-y)f(y)dy. (2.1) 

J — OO 

We can now follow the method of Walsh [18] to see that that (1.2) admits a 
mild solution u if and only if u is a predictable process that satisfies 



/CO ft 
/ b(u(s, y))p(t, y - x) ds 
-oo JO 

/'OO ft 
/ cr(u(s, y))p(t, y-x) w(ds dy). 
-oo J0 



dy 



(2.2) 



Here the proof of local existence follows via standard localisation means. The 
proof of our theorem will involve some bounds on the transition density function 
appearing in (2.1). Our first proposition recalls a few well known properties 
which we will use to derive those necessary bounds. We do not present proofs 
but mention [15], where the reader can find lots of details about symmetric 
stable processes. 

Proposition 2.1. The transition density p(t, ■) of a strictly a-stable process, 
satisfies the following 

(a) 

p(st, x) = r^ a p(s, t^^x), (2.3) 
(b) For t large enough such that p(t, 0) < 1 and a > 2, we have 

p(t, (x - y)/a) > p(t, x)p(t, x) for all x e R, (2.4) 

(c) 

p(t, x) x t-V- A (2.5) 



4 



As a first application of the above proposition, we have the following lemma. 
We will use this to derive a very useful ODE. 

Lemma 2.2. Let < s < c\t, where c\ is a strictly positive constant. Then the 
following inequality holds 



/ p 2 (s, x)p 2 (t, x)dx > c 2 t~ 
JR 



3/a 



where c 2 is a positive constant depending on a and c\ . 

Proof. Holder's inequality with respect to the probability measure p(s, x)dx 
gives 



J p(t, x)p(s, x)dx < p 2 (t, x)p(s, x)dx 



1/2 



Another application of Holder's inequality yields 

J p 2 (t, x)p(s, x)dx < p 2 (t, x)p 2 (s, x)dx^j p 2 (t, x)dx^j 

= \\p(t, -)IIl»(r) (^p 2 ^' x )p 2 ( s > x ) dx ^j > 

where \\p(t, •) lli 2 (R.) := JnP 2 (t> x)dx. Combining the above two inequalities, 
we obtain 



L 



(j R p(t, x)p(s, x)dx) 4 



p 2 (t, x)p 2 (s, x)dx > |2 



Note that 



/ p(t, x)p(s, x)dx = p(t + s, 0), 



and 

lb(*,-)||i 2( R)=p(2t, 0). 

Using the fact that p(t, 0) x i -1 /", we obtain the required inequality. □ 

Our next result concerns the blow-up of solutions to the class of ODE men- 
tioned above. We will use this to conclude the proof of our main result. 

Proposition 2.3. Consider the following ordinary differential equation. 

y'(s)s b = y(s) 1+a 

where a, b > with initial condition y(to). If b < 1, the solution blows up for 
any positive initial condition y(io). However, if b > 1, then the solution blows 
up whenever y(t ) is large enough. 
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Proof. We begin by looking at the simplest case when 6=1. Elementary 
calculus show that the solution y(t) satisfies 

y(t)- a = y(t )- a -aHl). 

to 

The above shows that y(t) blows up in finite time for any positive initial con- 
dition y (to). We now turn to the case when 6^1. Again some calculus shows 
that y(t) satisfies 

y(tr a = y(to)- a + ^^^ . 

In the case that b > 1, the above implies that y(t) blows up only when 

y(t y > (2.6) 

When b < 1, the solution blows up for any positive initial condition. □ 

We end this section by recalling Jensen's inequality which we will use fre- 
quently in the forthcoming section. 

Lemma 2.4. Let p(dx) be a probability measure on R and u be a non-negative 
function. Then the following holds whenever f is a convex function, 

f [J u(x)p(dx)) < Jj{u{x))p(dx). 

3 Proof of the main result. 

The proof of the main theorem hinges on the next proposition. We first fix a 
couple of notations. 
Denote 

F(t) := / E\u(t, y)\ 2 p 2 (t, y)dy, (3.1) 

JR 

and 

G(t) := f E\u(t, y)\p(t, y)dy, (3.2) 
Jr. 

where {u(t, x)} t >o, xen is the solution to (1.2). We can then state the following 
proposition. Sugitani [17] was the main source of inspiration for its proof. 

Proposition 3.1. The following two statements hold. 

(a) F(t) blows up in finite time if and only if for all ieR, E\u(t, x)\ 2 blows 
up in finite time, 

(b) G(t) blows up in finite time if and only if for all x 6 R, E|w(i, x)\ blows 
up in finite time. 



G 



Proof. We will only prove part (a) of the proposition. The proof of the second 
part is similar. Note that it is enough to show that E\u(t, x)\ 2 — oo whenever 
F(t) = oo for some large t. The other implication is easy. Using property (2.3), 
we can write 

Let c a := 2^+1- We can now use (2-4) to bound the right hand side of the above 
display as follows. 

S 



p(t- s, x -y) > \J— S ) P( s > x M s >y)> ( 3 - 3 ) 

for t < s < c a t 7 where t is chosen large enough so that p(s, 0) < 1 for all 
s > to. From the mild formulation of the solution i.e (2.2) and the positivity of 
b and the initial condition uo{x), we obtain 

E|u(t, x))\ 2 > f [ E\a(u(s, y))\ 2 p 2 {t-s, x-y)dyds. (3.4) 

JO JR 

Since the integrand appearing in the right hand side of the above display is 
strictly positive, we have 

E|«(t, x))\ 2 > I" f E\a(u(s,y))\ 2 p 2 {t-s,x-y)dyds. 
Jt Jr. 

Using the bound (3.3), and the lower bound on a together with Jensen's 
inequality, we can write 

E|«(t, x))\ 2 >^ " (j^) 1 P 2 ( s ' x ) J r VWHs, y))\V(s,y)dyds 

- Cl [ at ' p 2 (^)/ r ( e K^)I 2 ) 1+ V( S , y)dyds 

r=* ( s \ 2/a 

~ Cl J t {t~sj P^^ X )WP^ -)H 2 L 2 (R) 
x / (E\u(s, y)\ 2 r? J^l ? dyds 

rc a t / x 2/q 

- Cl i {— 8 ) p 2 ( s ^)\\p^-)Whw 



»(R) 

1+/3 



2 p 2 ( s > y) 

Lit/ J us. 

R) / 

(3.5) 



The result now follows from the above inequality and the strict positivity of the 
transition function. □ 
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We are now ready to prove our main result, Theorem 1.3. 

Proof. We proceed in two steps. The first step will be carried out assuming 
condition 1.1, while the second step will be completed under condition 1.2. 

Step 1: In view of Proposition 3.1, we aim to show that F(s) blows up in 
finite time. We start with (2.2). Bearing in mind that b is positive, we end up 
with 



/ E\u(t,x)\ 2 p 2 (t, x)dx > l 

JR JR 



u (y)p(t, x-y)dy 

R 

t 



2 

2/ 



p (t, x)dx 

III E\a(u(s, y))\ 2 p 2 (t - s, x - y)dsdyp 2 (t 1 x)dx 

JR JR JO 
■=h+h- 

We bound the I\ first as this is straightforward. The lower bound on the initial 
condition yields 

/ Uo(y)p(t, x - y)dy > n. 

JR 

Hence, we have 

h>K 2 \\p(t,-)\\h {R y 

To bound li , we will use some of the estimates from the proof of Proposition 
3.1. We first use inequality (3.5) to obtain 

/ E\u(s,y)\ 2 f^ y J dy) ds. 

jr \\p( s ' •)II L 2 (r) J 

We now use the fact that \\p(s, -)IIl 2 (r) x s_1 ^" an d the lower bound obtained 
in Lemma 2.2 to write 



h > cijT" (j^j t- 3 / a s^ a ^E\u(s,y)\ 2 p 2 (s,y)dy^ ds 

rc a t (2+/3)/a / r \ 1+/3 

> ci jf t5/a ^E\u(s,y)\ 2 p 2 (s,y)dyj ds 



We now combine the estimates for ii and I2 to obtain 



/ E\u(s,y)\ 2 p 2 (s,y)dy 

JR 

f c a t (2+/3)/a / r \^ 

" ° 2 It Ur E| " (S ' y))l V(s ' y)Ay ) AS 
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The definition of F and some algebra yield 



F{t)t 5 ' a > K 2 c 3 t 4 / a + c 2 s^+^/ a F{ S ) 1+f3 d S , . 
We set Y(t):= F(t)t b ' a to obtain 

1-Cc.t y( \i+p 

Y(t)> K *c 3 t^ + c 2 J^ds. 

J to 



We are interested in the blow-up of Y(-), by a comparison principle; see for 
instance [16], it suffices to consider the following ordinary differential equations. 

Y'(t)t<- 3+4 M a = C4Y(t) 1+f3 for t>t , 

with initial condition Y(t a ) = c 3 K 2 t^ a . Note that 3 + 4/3 > a, so we can 
use Proposition 2.3 to conclude that the above equation blows up in final time 
whenever k is large enough. Hence F blows up in finite time as well. An 
application of Proposition 3.1 concludes the first part of the proof. 

Step 2: As in the previous step, we start with the mild formulation of the 
solution but take expectation to obtain 



Eu(t,x)= / u (y)p(t, x-y)dy + / / Eb(u(s, y)p(t - s, x - y) dy ds 

We aim to show that the function G(s) defined by (3.2) blows up in finite 
time. Proposition 3.1 and an application of Holder's inequality will then yield 
the result. We begin by multiplying the above display by p(t, x) and integrate 
over R to obtain 

G(s) = u (y)p(t, x - y)dyp(t, x)dx 

Jtl Jtl 

+ / / / Eb(u(s,y))p(t-s,x-y)dydsp(t,x)dx 
JtlJo Jtl 

■■= h+h- 

I 3 is easily seen to be bounded below by k. We now bound I4. As in Proposition 
3.1, there exists a to > large enough so that for to < s < c a t, (3.3) holds. This 
bound together with the growth condition on b gives 



I4 = / / / Eb(u(s, y))p(2t — s, x — y) dydsp(t, x) dx 
Jtl J a Jtl 

> ^J^JJ (2^) J^E\<s,y)\ 1+ ^p{ S ,y)dydsp{t,x)dx. 
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We now apply Jensen's inequality twice to obtain 

cc a t / \ l/a / f \ 1+7 



h > c h JJ ^nu(s,y)\p(s,y)dy^ ^ ds 

S "l (j) G(s) 



G( o>. + c,p' (G( t 1/ T^ 



Combining the bounds on 7 3 and I4, we obtain 

(G(s)s 1 / a ) 1 +^ 

'to 

which in turn yields 

G(t)t^ a > nt 1 /" + c 5 (G(S ^' / / 1 )1+7 d S for i>t . 

We set Y(t) := G(t)t^ a and look at the following ODE; Y'{t)V' a = c 6 Y(t) 1+ ^ 
with Y(t ) = nty a . Now if 7/a < 1, then F(t) blows up in hnite time, when- 
ever the initial condition satisfies (1.3). But if j/a > 1, then Y(t) blows up 
whenever k is large enough. Since G{t) blows up whenever Y(t) blows up, the 
second part of the theorem is proved. □ 

We end this section with the following corollary whose proof is obtained by 
setting b{u) = in Theorem 1.3. 

Corollary 3.2. Consider the following equation 

d t u(t, x) = -(-A) a / 2 u(t, x) + a(u(t, x))w(t, x) 
u(0 , x) = u (x), 

Suppose that condition (1.1) is satisfied, then the second moment of the 
solution to (3.6) blows up in finite time as long as k is large enough. 

4 Remarks and Extensions. 

In [17], the author considered the following deterministic equation, 

d t u(t, x) - -{-A) a ' 2 u{t, x) + f(u(t, x)) 
u(0,x) = u (x), 

and explored non-existence of global solutions under certain suitable conditions 
on / and the initial data Uo(x). More precisely, it was shown that if f(x) > 
c\x\ 1+ P with < (3 < a, then the solution to the above equation blows up 
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for any positive initial condition. We can randomise the above equation by 
multiplying the forcing term by a white noise and hence end up with (1.2) 
but with b(u) = 0. We can then make the following interesting observation. 
According to Corollary 3.2 finite time blow up occurs for large enough initial 
data even for f3 < a. However, in the deterministic shown in [17], if 

(3 < a, there is blow up for any initial data. Thus in this case, one might say that 
the presence of noise might hinder blow up. That is there is a range of initial 
data, for which there is definite blow up in the deterministic case, but there 
might not be blow up in the stochastic case. It is also interesting to note that 
that the same mechanism has been reported in the literature for the stochastic 
wave equation, when forced by a Wiener field [10]. 

We end this paper with a couple of possible extensions. Our main result 
was proved under the assumption that d = 1. It would be interesting to prove a 
similar result for d > 1 but with a noise which has some spatial covariance. One 
would thus need to find suitable conditions on the spatial covariance as well. 
This might not be an easy task as has been shown in [8] where such equations 
were studied. One should also note that under Condition 1.2, blow up mainly 
occurs due to the additive nondinearity and the theorem is proved by looking 
at the expectation of the solution. Hence in this case, one can easily extend the 
result to higher dimensions. 

Another interesting question is to consider blow up or global existence of 
(1.2), when b(u) plays the role of a damping term, instead of a source term. In 
this case there will be competition between b(u) and a(u)w(t, x), and possible 
blow up will ultimately depend on which term dominates. 
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